Stability of Cosmological Solution in f{R) Models of Gravity 
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We reconcile seemingly conflicting statements in the literature about the behavior of cosmological 
solutions in modified theories of gravity where the Einstein-Hilbert Lagrangian for gravity is modified 
by the addition of a function of the Ricci scalar, f{R). Using the example of f{R) — ^ijf^ /R we 
show that only such choices of f{R) where f /AP? > have stable high-curvature limits and well- 
behaved cosmological solutions with a proper era of matter domination. The remaining models enter 
a phase dominated by both matter and scalar kinetic energy where the scalar curvature remains 
low. 



The addition of a new term f{R) to the Einstein- 
Hilbert action provides a possible explanation for the 
late-time acceleration of the universe [H, 0, However 
there is some confusion in the literature as to what kind 
of functions f{R) result in solutions that are stable and 
cosmologically viable. In Q, the authors claim that a 
cosmology which behaves like ACDM at early times and 
is stable to perturbations requires //jfl, = d^f/dR^ ^ 0. 
On the other hand, Bean et al. [sj have demonstrated 
that linear perturbations are stable in the CDTT model 
[H for which < 0. In addition, Amendola et al. @, 0| 
have claimed that f{R) models do not have a normal 
matter-dominated epoch. 

The aim of this Brief Report is to clarify the situation. 
We show that the sign of /rr does determine whether 
the theory approaches the GR limit at high curvatures. 
For Jrr^ > 0, the theory at high curvatures behaves very 
close to GR and is stable. For Jrr < 0, the GR limit 
is unstable. In an f{R) model, high density does not 
necessarily correspond to high curvature. We will show 
that the Ricci scalar rapidly evolves to a low curvature 
solution at high redshift in this case. This is the insta- 
bility discovered in [3|. However, the low-curvature solu- 
tion itself is stable, and hence it has well-behaved linear 
perturbations The background cosmology in this low- 
curvature phase does not have the usual period of matter 
domination with H oc oT^I"^ . Instead it enters a (/)MDE 
phase @, 0| where H oc a~^, and is not observationally 
viable. 

Preliminaries — We consider modifications to the 
Einstein-Hilbert action of the form Isl 
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where Rj is the Ricci scalar, = SirG and Cm is 
the matter Lagrangian. We define Jr = d//di?,j and 
/rr = d^f/dRj. "J" denotes the Jordan frame, to be 



distinguished from the Einstein frame below. We con- 
sider only the metric-formalism version of the theory, in 
which the Ricci scalar is a fully dynamical degree of free- 
dom. 

In order for an model to be consistent with cos- 

mological observations, the evolution of the scale factor 
oj must remain close to that of the concordance ACDM 
model: i.e. a period of radiation domination, followed by 
high-curvature matter domination, ending in an era of 
acceleration. In Q, the authors showed that there exists 
an f{Rj) that reproduces any desired expansion history. 
The function is parameterized by 
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where Hj is the Hubble parameter. However, models 
with B < were found to be unstable to linear perturba- 
tions at high curvature and hence to not reproduce the 
desired phenomenology. 

We compare this finding with that in the Einstein 
frame in which stable cosmological solutions for B < 
models have been found 0. As emphasized by in 
the Einstein frame gravity is unmodified but there is an 
additional minimally coupled scalar which evolves in a 
potential. Provided that /rr ^ and H-Jij > 0, we can 
use the conformal transformation to the Einstein metric, 
gjf^ = exp(o-0)5uiy ■^ith the new scalar field defined as 



V; = fTln(l + /fl) 
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where a — sgn(/ij) ensures that V 5^ 0- The Jordan- 
frame curvature Rjiip) can then be found by inverting 
the relation The resulting Einstein-frame action is 
one for a scalar-tensor theory of gravity [3| 
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"E" denotes variables in the Einstein frame. Here the 
matter is not minimally coupled, which results in the 
need to transform back to the Jordan frame for the inter- 
pretation of physical results. Henceforth all derivatives, 
with the exception of Jr and /rr, will be taken with 
respect to Einstein-frame variables. 

High curvature stability — For expository purposes we 
use a generalization of the CDTT model [l|: 
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Assuming that tp <^ 1, i.e. that the curvature is high — 
i?j /i^ — we can expand the effective potential around 

V' = 0: 
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This potential has an extremum at i/'gr ~ M^/('^^Pm)^j 
i.e. at the GR limit of i?j « K'^pi^. For the CDTT model 
this is a maximum, whereas for the mCDTT model this 
is a minimum. The potentials for both models are pre- 
sented in Fig. [TJ 



where a = ±1 and /i ~ Hq. The a — +1 model is 
the original CDTT choice with B < and has a pe- 
riod of late-time acceleration with an effective w — —2/3. 
As in all inverse curvature models, the behavior in the 
low curvature acceleration phase where {/rI ^ 1 is in- 
dependent of a [ll[. However the acceleration phase for 
cr = — 1 (mCDTT) eventually causes a sign change in 
1 + /fl evolving past a coordinate singularity in the Ein- 
stein frame. Moreover, a sign change in 1-1-/^ and hence 
B at low curvature would lead to Jordan-frame instabil- 
ities in a perturbed universe. Nevertheless the mCDTT 
model is illustrative and its pathologies are avoidable. 
Any function f{R) such that limii;j_,oo fR = could have 
been chosen for this analysis. For example, the model 
f{R) = ^2(i?/^2)m ^-^^^^ < TO < 1 has 5 > and ex- 
hibits stable and observationally viable acceleration [l^ . 

We analyze the model in the Einstein frame, with the 
scalar defined as in At high curvature, ip is inversely 
related to the Jordan-frame Ricci scalar i?j 



lim tp 

-Rj— ^oo 



(R^y 



(7) 



Small tp implies a high curvature. If ip evolves away from 
zero, the expansion enters a low-curvature phase. The 
potential for the scalar, as defined in ^ is 



cr(e^'> -1) . 
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We assume that the matter is a sum of a pressureless 
component, Pm — 0, with density pm and an ultrarel- 
ativistic component, Pi- — 1/3, with density pr- Non- 
minimal coupling to the trace of the energy-momentum 
tensor requires that the evolution of ip is driven by an 
effective potential 
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where we have used p 
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g2o-i/' bring out the depen- 



dence on tp. Neither of the two quantities is conserved in 
the Robertson- Walker background of the Einstein frame. 
For example, p^ evolves according to 
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Figure 1: Exact effective potentials with constant term sub- 
tracted for a range of matter densities for scalar field in CDTT 
(top panel) and mCDTT (bottom panel) models. Both mod- 
els exhibit extrema at -Rj = K^p'^- a maximum for CDTT 
and a minimum for mCDTT. 



We can easily demonstrate the instability of the high- 
curvature solution for the CDTT model in the Einstein 
frame. For ^p 'ti 1, the equation of motion of the scalar 
field becomes ([9]) 



where ' 
tion 
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d/dlnoE, and obeys the Friedmann equa- 
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We now assume that there exists a well behaved high- 
curvature solution to (fT2|) . say ipo. We then perturb 
around this solution by defining tp — ipo + Sip which re- 
sults in the linearized equation of motion for the pertur- 
bation 
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For a = +1, Sip is unstable on very short timescales 
since i/'o 1- However, for the mCDTT model the high- 
curvature solution is stable, provided it exists. This anal- 
ysis is related to where the mass of the scalar around 
= 



=0 — /rr - 3i? - 4/ , 



(15) 



gave a large tachyonic value < for the CDTT model 
at high curvature. The field rapidly rolls away from this 
maximum to high values, causing the Jordan frame cur- 
vature to be significantly below that of GR. We illustrate 
this with a numerical solution in Fig. [21 The calculation 
is initialized at the GR value of the Ricci scalar but it 
drops quickly to a small value and proceeds to decay as 



10 



10 



10 



a. 



10^ 



10" 



10' 




-15 



Figure 2: Evolution of the Jordan-frame Ricci scalar for 
LCDM (solid; fi = Ho), mCDTT (dashed; a = 2) and CDTT 
(dotted; a = 2). All models are started with the same GR ini- 
tial conditions. This solution is unstable in the CDTT model 
and the curvature rapidly relaxes to the <^MDE phase where 
it is much lower than for the corresponding GR solution. This 
phase replaces the standard matter domination for f{R) mod- 
els with B < 0. mCDTT follows standard LCDM evolution 
until late times, where it asymptotes to a milder acceleration 
phase. 



We now need to estabHsh the existence of the high- 
curvature solution in the case of mCDTT, i.e. that cor- 
rections to GR remain small throughout the matter- 
dominated phase. We continue to assume that V', V'' ^ 1 
(implying that Jordan- and Einstein-frame quantities are 
approximately equal) and we ignore the contribution of 
radiation. We take tpGR — /^^/('^^Pm)^ the GR so- 
lution for the Jordan-frame Ricci scalar. Defining the 
perturbation away from this solution as 5ip = ip — ipGR, 
we obtain the linearized equation for these perturbations 



from (fT2)) during matter domination 



H = -45'(/'GR 
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a" = 'i^^ / plJ^aj — 1). This is a damped 
harmonic oscillator with an extremely high frequency, 
Lo ~ 1/a^ag. The amplitude is driven by the inhomo- 
geneous term, giving 
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This remains extremely small throughout matter domi- 
nation, confirming that for the mCDTT model the cos- 
mology remains extremely close to that for GR until the 
acceleration phase. 

Low-curvature (j)MDE phase — We have demonstrated 
above that in the B > mCDTT model the high- 
curvature solution is stable, whereas it is not in the B <Q 
CDTT model. The evolution of f{R) theories generically 
exhibits a fixed point in which the Einstein-frame scalar 
field is kinetic-energy dominated and remains a constant 
fraction of the energy density, together with the matter 
This phase replaces the ordinary matter dominated 
expansion for B < and is called the </)MDE solution. 
The analysis was extended in [l3| to show that B > Q 
models do have usual matter-dominated eras for at least 
some initial conditions. 

In this section we will demonstrate — using a slight 
modification of the method originally introduced in [l3| 
and then used in — that indeed the 0MDE phase is 
a saddle point in the evolution for both CDTT and 
mCDTT, but it can only be achieved for CDTT. 

We introduce the dimensionless variables 
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Assuming that radiation is negligible, we can rewrite the 
Friedman equation lfT3|) as an equation for the matter 
energy density and its derivative 
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We assume that the potential for the scalar field is of ex- 
ponential form, Viijj) = Aexp(7'0), which is valid when 
V' > 1. From lEl), 7 = -3/2 for CDTT, while 7 = 2 
for mCDTT. This allows us to recast ^ and ^ as 
evolution equations for x and y: 



2^ = 77 ^ 7:x — —X 
2 2 2 
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The fixed points of the evolution occur wherever x' = 
y' = 0. In particular, the point corresponding to 0MDE 
lies at {x,y) = (it/3,0). At this point, the scalar field 



4 



energy density is completely dominated by the kinetic 
term and the fraction of the total energy in matter is 
= 8/9. From equation ifTO]) . the effective equation of 
state is w = 1/9, resulting in the Einstein-frame Hubble 
parameter evolving as (x a-^^^^. The Jordan and 
Einstein frame Hubble parameter are related by 

= He 0-^^/2 (^1-^) . (22) 

Since ^ 21naE/3 at this fixed point, Hj oc aj^, just 
as during radiation domination. 

Using l|12p and (fT3| . we can express the Jordan-frame 
Ricci scalar as 

R,=e^^(iV + 2a^^ , (23) 

and hence both the models under consideration have 
i?j < at the (j)MBE fixed point. 

The eigenvalues at this point are —4/3 and 2/3(5+ao'). 
For both models 0MDE is a saddle point. However, only 
for the CDTT model is this stage of evolution possible to 
achieve. Here, the field has a large value and is running 
off to positive infinity, exponentially suppressing the con- 
tribution of the potential to the Hubble parameter. In 
the mCDTT model, V (xe'^'f': the large-field limit is also 
the large-potential limit. The existence of the 0MDE 
phase therefore requires that V' ^ 0, which is not a valid 
limit of the Jordan-frame theory and is a result of the 
use of the exponential approximation to the potential. 

Discussion — We have shown that in the case of the 
generalized CDTT model, the sign of B determines the 
behavior of the model and leads to a qualitative change 
in the Einstein-frame potential for the scalar field. This 
scalar field has a large negative mass squared in the high 
curvature limit for B < models, resulting in unstable 
solutions in the GR limit. The field runs off to a large 
value, where the potential is effectively exponential. The 
cosmology enters a phase of simultaneous scalar field and 



matter domination with the physical Hubble parameter, 
Hj, proportional to aj^ and a very small value of the 
curvature scalar in the Jordan frame. 

The situation is radically different in the case of 
mCDTT, where B > at high curvature. The GR-like 
solution is stable, the cosmology goes through a phase of 
standard matter domination where corrections to GR re- 
main small. The acceleration behavior is asymptotically 
the same as CDTT but this specific model will develop 
instabilities deep in the acceleration regime. 

Another important difference between the two model 
classes parameterized by B is their behavior around mass 
sources, i.e. their predictions for Solar-System tests. In- 
deed, in i? < models the Ricci scalar is perturbed only 
slightly from its low background value by sources of mat- 
ter density, Rj ^ n^Pm- PPN parameters can then be 
calculated using the weak-field approximations [H, 
and which results in a highly excluded value of 7 = 1/2 
[lo| . However, this is not necessarily true for B > mod- 
els. For these models, the value of the Ricci scalar can 
be high, refiecting the local matter density. It is there- 
fore inappropriate to linearize the perturbations to 
As shown in and the chameleon mechanism re- 
sulting from the non-minimal coupling of the theory to 
matter and the existence of a minimum in the scalar po- 
tential around the GR value of the Ricci allows for the 
theory to evade Solar-System tests, at least for certain 
functions f{R) and certain ranges of parameters. We 
will elaborate in detail on the constraints put on f{R) 
theories by solar-system tests of gravity in a future work 




Acknowledgments — We thank R. Bean, S. Carroll, H. 
Peiris, J. Santiago, Y.S. Song, J. Weller and A. Upadhye 
for useful conversations. This work was supported by 
the KICP under NSF PHY-0114422, DOE under DE- 
FG02-90ER-40560, and the David and Lucile Packard 
Foundation. 



[1] S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, 
Phys. Rev. D70, 043528 (2004), astro-ph/0306438. 

[2] S. Capozziello, S. Carloni, and A. Troisi (2003), astro- 
ph/0303041. 

[3] S. Nojiri and S. D. Odintsov, Phys. Rev. D68, 123512 

(2003), hep-th/0307288. 
[4] Y.-S. Song, W. Hu, and I. Sawicki, Phys. Rev. D75, 

044004 (2007), astro-ph/0610532. 
[5] R. Bean, D. Bernat, L. Pogosian, A. Silvestri, and 

M. Trodden, Phys. Rev. D75, 064020 (2007), astro- 

ph/0611321. 

[6] L. Amendola, D. Polarski, and S. Tsujikawa, Phys. Rev. 

Lett. 98, 131302 (2007), astro-ph/0603703. 
[7] L. Amendola, D. Polarski, and S. Tsujikawa (2006), astro- 

ph/0605384. 



[8] A. A. Starobinsky, Phys. Lett. B91, 99 (1980). 
[9] G. Magnano and L. M. Sokolowski, Phys. Rev. D50, 5039 
(1994), gr-qc/9312008. 

[10] T. Chiba, Phys. Lett. B575, 1 (2003), astro-ph/0307338. 

[11] O. Mena, J. Santiago, and J. Weller, Phys. Rev. Lett. 
96, 041103 (2006), astro-ph/0510453. 

[12] T. Faulkner, M. Tegmark, E. F. Bunn, and Y. Mao 
(2006), astro-ph/0612569. 

[13] L. Amendola, R. Gannouji, D. Polarski, and S. Tsu- 
jikawa, Phys. Rev. D75, 083504 (2007), gr-qc/0612180. 

[14] E. J. Copeland, A. R. Liddle, and D. Wands, Phys. Rev. 
D57, 4686 (1998), gr-qc/9711068. 

[15] A. L. Erickcek, T. L. Smith, and M. Kamionkowski, Phys. 
Rev. D74, 121501 (2006), astro-ph/0610483. 

[16] T. Chiba, T. L. Smith, and A. L. Erickcek (2006), astro- 



5 



ph/0611867. [19] Defining <p = crip / ^/2/2>k yields canonical normalization. 

[17] I. Navarro and K. Van Acoleyen (2006), gr-qc/0611127. 
[18] W. Hu and I. Sawicki, Phys. Rev. D submitted (2007), 

arXiv:0705.1158 [astro-ph]. 



